From single particle to superfluid excitations in a dissipative polariton gas 
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Using angle-resolved heterodyne four- wave- mixing technique, we probe the low momentum ex- 
citation spectrum of a coherent polariton gas. The experimental results are well captured by the 
Bogoliubov transformation which describes the transition from single particle excitations of a normal 
fluid to sound-wave- like excitations of a superfluid. In a dense coherent polariton gas, we find all the 
characteristics of a Bogoliubov transformation, i.e. the positive and negative energy branch with 
respect to the polariton gas energy at rest, sound-wave-like shapes for the excitations dispersion, 
intensity and linewidth ratio between the two branches in agreement with the theory. The influence 
of the non-equilibrium character of the polariton gas is shown by a careful analysis of its dispersion. 



PACS numbers: 67.10.-j, 71.36.+C, 78.47.nj, 78.67.De 

The demonstration of superfluidity of Helium II by 
Kapitsa [1 , Allen and Misener [2 in 1938 triggered fun- 
damental theoretical research for understanding quan- 
tum fluids. London quickly linked superfluidity with 
Bose Einstein condensation [3 , stressing the importance 
of the bosonic nature of the particles in the phenomenon. 
Meanwhile, Landau developed the idea of sound-like ex- 
citations of a superfluid [4]. These intuitions were later 
confirmed by Bogoliubov whose microscopic theory of the 
weakly interacting Bose gas [5 reveals a superfluid phase 
in which elementary excitation is a coherent superposi- 
tion of counter-propagating particle and hole. 
In a superfluid at equilibrium, the momentum dispersion 
of the excitations should deviate substantially from the 
parabolic single particle case. In addition to a global 
blueshift due to repulsive particle-particle interactions, 
the dispersion turns linear for small wavevectors as a re- 
sult of the sound-wave nature of the superfluid elemen- 
tary excitations. In addition to the normal positive en- 
ergy branch (NB), the Bogoliubov dispersion is expected 
to display a negative energy "ghost" branch (GB), result- 
ing from the hole component of the excitation [6 . The 
NB and GB are the mirror image of each other and rep- 
resent the dispersion of Bogoliubov excitations [7 . Su- 
perfluid excitations have been extensively investigated in 
both He superfluids [SHTO] and ultracold atomic gases 
pTHT3] . First signs of a Bogoliubov transformation have 
been given by Vogels et al [14 in atom condensates. 
In a semiconductor microcavity where quantum well ex- 
citons and cavity photons are strongly coupled, light- 
matter quasi-particles called polaritons are formed. 
These composite bosons show unique properties which 
promote them as a model system of the weakly inter- 
acting Bose gas. Indeed, since demonstration of Bose- 
Einstein condensation of polaritons [15], superfluidity of 
coherent polariton fluids has been assessed the observa- 
tion of zero-viscosity and critical velocity by Amo et al 
[16l [IT] . Very recently superfluid excitations have been 



investigated in a nonresonant polariton condensate [TE\ . 
The negative Bogoliubov dispersion branch could surpris- 
ingly not be observed though. 

However, the peculiarity of polaritons with respect to 
other superfluids is their non equilibrium nature. In the 
case of non-resonant excitation of the polariton conden- 
sate the theory predicts the appearance of a flat diffusive 
region in the Bogoliubov dispersion around k=0, result- 
ing from the coupling between the exciton reservoir and 
the polariton field [7 . Under cw resonant excitation no 
exciton reservoir is created and the excitation spectrum is 
therefore much closer to the equilibrium Bogoliubov spec- 
trum [19 . Under pulsed resonant excitation however, we 
expect strong modifications in the polariton fluid disper- 
sion due to the decay of the polariton population. 
In the present work, following the suggestion of a recent 
theoretical work [20] , we choose the following strategy to 
investigate the superfluid character of a dissipative coher- 
ent polariton gas (CPG). We use the four- wave mixing 
(FWM) technique to stimulate emission from both the 
normal and negative energy branch in the dispersion of 
a resonantly created CPG (Figjl^). Our experimental 
access to both these dispersions allows us to observe the 
modification of the polariton gas excitations and to con- 
clude on the superfluid character of the CPG, e.g. Bo- 
goliubov transformation. Contrary to reference [20 we 
use circularly polarized pulsed excitation and thus probe 
the excitations of a spin polarized CPG in the dissipa- 
tive regime. We can analyze the dynamics of the CPG 
excitations as the polariton density decays. We evidence 
the influence of dissipation on the excitation spectrum of 
a CPG resulting in an asymmetry of the NB and GB in 
the dispersion and an overestimated sound velocity. 
We investigated a single 8 nm Ino^o^Gao^geAs quantum 
well embedded in a GaAs wedge-shaped-cavity [21] . The 
Rabi-splitting is 3.4 meV at zero detuning. The sample 
is kept in a continuous flow cryostat at around 5K. We 
use a two pulse heterodyne FWM setup. For optical ex- 
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citation a 80 MHz pulsed Ti: Sapphire laser is used. The 
250 fs pulses are centered on the lower polariton k=0 en- 
ergy and cover the energy range of both lower and upper 
polariton dispersion curves. The laser beam is split in 
three parts: a reference, a pump(index 1) and a trigger 
(index 2). These last two co-circularly polarized beams 
are radio-frequency shifted with acousto-optic modula- 
tors (AOMs) by uji (75 MHz) and UJ2 (79 MHz) and are 
focused on a 100/im diameter spot on the sample. The 
pump pulse generates a CPG at rest (/ci=0). In a second 
time the trigger pulse stimulates excitations in the CPG. 
It is ten times weaker than the pump pulse and carries 
a finite wavevector k, which is varied between A:2 = to 
1 iim~^. Stimulated parametric scattering of polaritons 
from the CPG results in the generation of a FWM signal 
at opposite wavevector — /c2. The FWM signal is selected 
at kFWM= 2/ci — k2=—k2 and is directed into a mixing 
AOM together with the reference. The AOM, driven at 
ujfwm='^^i — ^2 (71MHz), up-shifts the diffracted ref- 
erence field and down-shifts the diffracted signal field by 
^FWM- This results in a spectral overlap of reference and 
signal field and a 7r-phase shift for the two heterodyne 
channels. The mixed beams are then dispersed in a spec- 
trometer. To recover the interference term between the 
FWM signal and the reference field, we subtract the two 
TT-shifted interferograms. In this way the classical noise, 
e.g. photoluminescence, is largely suppressed. From the 
measured interference we retrieve then the signal in am- 
plitude and phase by spectral interferometry [22|. 
Using nonlinear spectroscopy off-resonance polariton 
scattering has already been observed by Savvidis et al. 
|23j . It resulted in multi- wave- mixing processes between 
pump, signal and idler modes of a polariton OPO, reveal- 
ing non negligible nonlinearity in the microcavity. How- 
ever, the experiment could not allow them to demon- 
strate any modification of the dispersion in the vicinity 
of the CPG momentum. In our case, thanks to the het- 
erodyne detection scheme, we precisely obtain an exper- 
imental access to the low momentum excitations of the 
polariton fiuid. Heterodyning allows us to separate the 
FWM signal from linear coherent emission and higher 
order nonlinear emission of the CPG when the standard 
angular selection becomes inefficient. 

Bogoliubov transformation should result in a drastic 
change of the FWM response in terms of intensity con- 
tributions of the NB and GB modes as well as in terms of 
dispersion curvature. We therefore study the changes of 
the FWM spectrum with excitation intensity, and show 
that the nature of the CPG excitations changes from sin- 
gle particles to sound waves when increasing the polari- 
ton density. Any nonlinearity in the microcavity does 
result in the observation of a FWM signal on the nega- 
tive energy side. Here we show that this "off resonant" 
emission changes drastically when the superfluid regime 
is reached. In this case the emission intensities of NB 
and GB become comparable. 
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FIG. 1: a With a fs-pump pulse the coherent polariton gas is 
created. The trigger pulse gives rise to two distinct parametric 
scatterings (dashed or plain red arrows), b FWM spectrum 
(linear scale) as a function of the excitation photon density 
for a trigger wavevector k=l iim~^ . The white dashed line 
shows the k=0 polariton energy and the two white solid lines 
guide the eye and display the evolution of the NB and GB. 
A threshold at Pth— 5.10^*^ photons. pulse" ^cm~^ is observed 
indicated by the red dotted line, c, d and e show FWM 
spectra in the perturbative regime, short above threshold and 
the superfluid regimes respectively. The red curves display the 
Lorentz fits for the two resonances. 



In Figure [TJd, we plot the measured FWM spectrum as 
a function of excitation density at a fixed pump-trigger 
delay of 6.5ps. The trigger incidence angle was kept 
constant to stimulate excitations at wavevector 1 /im~^ 
where the NB and GB are well separated. Two emission 
lines are observed. One lays below the k=0 polariton en- 
ergy, the other one, laying above the k=0 energy, is the 
normal branch. At low excitation intensities (below the 
threshold Pth) the emission from the NB is much stronger 
than the emission of the negative energy branch (Fig. 
[l]^). Besides the emission of the negative energy branch 
is symmetric to the NB emission with respect to k=0. 
These two features demonstrate that the system is still in 
the perturbative regime, where the dispersion is given by 
the parabolic single-particle one. The emission at nega- 
tive energies is attributed to off-resonance polariton scat- 
tering processes. Experimentally we observe a threshold 
Pth for the GB intensity at 5.10^^ photons. pulse~^cm~^. 
Above Pth (dashed line) the GB emission is steeply en- 
hanced and becomes of the same order of magnitude as 
the NB emission. This shows that the excitations of the 
CPG have changed. Fig. [T]i shows the FWM spectrum 
just above Pth- Far above Pth (Pigfl^)^ the average in- 
tensity ratio between GB and NB is around 0.9. 
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The linewidth is also affected by the increase of the po- 
lariton density. In the low density regime the linewidth 
of the negative energy mode is much broader than the 
one of the NB mode. The Lorentz fits of the two reso- 
nances in FigjT]: (in red) exhibit a FWHM ratio between 
the GB and NB of 2.6 (FWHM: 0.14 meV(0.36 meV) for 
the NB (GB)). The change of the regime is accompanied 
with a redistribution of the linewidth. As its density in- 
creases and the ratio between the GB and NB intensity 
tends to one, the linewidths of the GB and NB become 
equal. Just above Pth (Fig [l]i ) the FWHM ratio has 
already decreased to 1.8 (FWHM: 0.17 meV(0.30 meV) 
for the NB (GB)). At high density (Fig. ^) the NB and 
GB linewidths are 0.26 meV and 0.27 meV, respectively. 
The relative intensities of the positive and negative en- 
ergy contributions of the FWM signal can be calculated 
in the lower polariton (LP) subspace, in which the po- 
lariton dynamics are described by the time-dependent 
Schrodinger equation: 



i^(r,t) = (e(V)+^|^(r,t)r 



z|)^(r,t) 



where ^ is the LP wavefunction, e its kinetic enegy, g 
the polariton interaction constant, and Fp^t the pump and 
trigger excitations. According to the design of our exper- 
iment, three modes must be considered. The CPG mode 
at k=0 is described by cpo. In analogy with Bogoliubov's 
theory, we call u and v* the counter-propagating per- 
turbative modes coupled by polariton interaction. The 
trigger pulse excites the v* mode and FWM is generated 
on the u mode. The LP wave function reads as: 



^(r, t) = (^o(r, t) + u{r, t)e'^' + ^*(r, t)e 



-ikr 



(2) 



Equation [T] can be solved numerically and the spectrum 
of the u mode can be calculated to obtain the intensity 
ratio between GB and NB. Fig|2^ provides the calculated 
density evolution for the emission contribution of NB and 
GB as a function of the normalized blueshift (gno/j) for 
different excitation wavevector. At low densities the cal- 
culated ratio is 1/3 for any wavevector. As the density in- 
creases, the elementary excitations are modified and the 
u and V modes couple and redistribute polaritons in both 
positive and negative energy branch. The ratio tends 
faster to one for smaller k as expected by the equilibrium 
Bogoliubov theory, because for smaller k the transition 
is more Bogoliubov like. The theoretical intensity ratio 
of 1/3 in the low density regime comes from different 
damping times of the FWM resonances: it is three times 
larger for the GB resonance than for the NB resonance. 
This is also what we can extract from the linewidth of 
the emission lines as shown in Fig. [T]^. These changes in 
the FWM spectrum demonstrate the appearance of the 
Bogoliubov excitations in the CPG. 
The excitation with fs-laser pulses generates populations 
of both LPs and upper polaritons (UPs). This yields os- 
cillations in the excitation intensity spectrum for the NB 
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FIG. 2: a Calculated intensity evolution of the NB (blue) and 
GB (red) emission contribution plotted for different pertur- 
bation wavevectors. Calculation parameters: h'j = O.lmeV, 
h/2m — 0.bmeV./j.m^ . b Normalized intensity of the NB and 
GB (blue and red crosses respectively and dashed lines for 
the fits). The blue and red squares display the normalized 
intensities of the NB and GB with removed oscillating part. 



as it can be seen in Fig[T]3. These oscillations complicate 
the discussion of the NB/GB intensity ratio. They could 
in principle be avoided by shaping the excitation pulses 
but this would not allow us reach the high intensities re- 
quired for this experiment. The normalized intensities 
of the NB and GB are presented in Figj2]3 as a function 
of excitation intensity. The oscillation imposed by the 
fs-pulses is clearly visible (crosses in ^p)- In order to 
suppress this feature due to the LP-UP mixing, we fitted 
the data with an combination of an exponentially decay- 
ing sine added to a standard exponential variation (dased 
lines in[2]3). The oscillating part is then subtracted from 
the data to obtain the mean GB (NB) intensity contri- 
bution (squares in|2]3). At low excitation power the NB 
emission is 4 times more intense than the GB emission. 
This is not far from the value of 3 expected by the pertu- 
bative theory. With increasing excitation power, the NB 
contribution decreases whereas the GB contribution in- 
creases. At high excitation power, the intensities of the 
two branches tend to the same value as it is expected 
for the Bogoliubov transformation. The calculated den- 
sity evolution for the NB and GB emission contribution 
(Figj2^) is displaying the experimental behaviour. 
To measure directly the curvatures of the NB and GB 
dispersion, we analyze the dependence of the FWM spec- 
trum on the trigger wavevector. Figure [3^ displays 
the spectra of the FWM signals obtained for different 
wavevectors of the trigger pulse integrated between 5 
and 6 ps delay. The intensity of the two pulses and 
quite importantly of the pump pulse (12.5 mW, 10^^ 
photons/pulse/cm^) was kept constant ensuring the same 
polariton density in the superfluid phase for all measure- 
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FIG. 3: a FWM spectra at different excitation wavevectors. 
b Energies of these lines (red dots) versus wavevector. We 
compare them with the single-particle parabolic dispersion 
(red-dashed line) and with the FWM signal calculated (black) 
using the Gross-Pitaevskii equations (19]. 



with the GB. The linearization of the Gross-Pitaevskii 
equation in the case of a decaying superfluid, assuming 
that the variations of the energies are slow, confirms the 
later result: the steeper slope of the NB originates from 
the time integration of the decaying polariton emission 
and not from a physical change of the sound velocity. 
Finally, one observes replicas of the NB at lower energy 
in the calculation which are also a result of the decaying 
superfluid density. 

In conclusion we have evidenced the Bogoliubov trans- 
formation at the transition from a single particle gas to 
a superfluid in a semiconductor microcavity. We demon- 
strated the modification of the polariton dispersion cur- 
vature from parabolic in the pertubative regime to the 
peculiar linear dispersion of a dissipative superfluid. De- 
spite their dissipative character, microcavity polariton 
gases can be qualitatively described in the frame of Bo- 
goliubov theory of superfluids, closely linking polariton 
phenomenology to cold atom physics. 

We acknowledge R. Houdre and U. Oesterle for pro- 
viding the sample and SNSF for financial support. 



ments. The energies of the two emission lines are plotted 
in Figure |3]3 versus the trigger wavevector. The parabolic 
dispersion of polaritons obtained from transmission mea- 
surements in the low-density regime is also displayed for 
comparison (red dashed line in ^jp) . Theoretically, in a 
steady state regime the dispersion of both the NB and GB 
is expected to be symmetric and linear ^ il9^ . Indeed the 
shape of the NB in our experiment is linear and strongly 
differs from the standard single-particle parabolic disper- 
sion. This is the consequence of polariton-polariton inter- 
actions and a clear indication of superfluidity. However, 
there is no symmetry between the GB and the NB at 
low k. Furthermore, if we extract a linear slope from the 
NB dispersion we find a sound velocity Cs = 0.6 fim/ps. 
It does not correspond to the velocity obtained from the 
k=0 blueshift of 0.15 meV with the equilibrium formula: 
Cs = \fgnjM. The apparent sound velocity is thus larger 
than expected. This peculiar dispersion originates from 
the dissipative nature of polaritons. Taking the Bogoli- 
ubov equations, two physical processes can be pointed 
out: a blueshift of the whole dispersion corresponding to 
the blueshift of the condensate itself and the coupling of 
the two counter propagating modes of the Bogoliubov ex- 
citations leading to the change of curvature. If the former 
is linear with the particle density the latter is sublinear. 
The change of Bogoliubov-like curvature is therefore less 
sensitive to the polariton decay than dynamical blueshift. 
This leads to a steeper slope of the NB resulting in an ap- 
parent higher speed of sound and to the observed asym- 
metry in the time- integrated results. A more complete 
theoretical study is provided in additional information. 
Our calculation reproduces properly the different features 
of the experimental dispersion: the linearity of the NB 
with an overestimated sound velocity and the asymmetry 
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